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ABSTRACT 

An investigation was undertaken to measure the absorption coefficient, 
at a wavelength of 10.6 Micrometers, for mixtures of water vapor and a 
diluent gas at high temperature and pressure. The experimental concept was 
to create the desired conditions of temperature and pressure in a laser 
absorption wave, similar to that which would be created in a laser propulsion 
system. A simplified numerical model was developed to predict the 
characteristics of the absorption wave and to estimate the laser intensity 
threshold for initiation. A non-intrusive method for temperature 
measurement utilizing optical laser-beam deflection (OLD) and optical spark 
breakdown produced by an excimer laser, was thoroughly investigated and 
found suitable for the non-equilibrium conditions expected in the wave. • 
Experiments were performed to verify the temperature measurement technique, 
to screen possible materials for surface initiation of the laser absorption 
wave and to attempt to initiate an absorption wave using the 1.5 kW carbon 

dioxide laser at The University of Tennessee Space Institute (UTSI). The OLD 

technique was proven for air and for argon, but spark breakdown could not be 
produced in helium. It was not possible to initiate a laser absorption wave 

in mixtures of water and helium or water and argon using the 1.5 kW laser at 

UTSI, a result which was consistent with the model predictions. 
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I. INTRODUCTION 


Laser energy beamed from a remote laser to an orbital transfer vehicle can pro- 
vide a propulsion system with an exceptionally high specific impulse [1,2]. Several 
systems studies [3,4,5] have identified missions for which a beamed laser propulsion 
system can have significant cost or performance advantages. Several different meth- 
ods have been proposed to absorb the laser energy [6,7], but molecular absorption 
and plasma absorption appear to be the most promising. 

An ideal absorbing medium is one which can absorb the entire laser beam 
energy in a relatively small volume, and then be expelled from the vehicle as the 
propellent. In addition to these essential requirements, the absorber should have 
a low molecular weight to provide high specific impulse and be relatively easy to 
store on-board the vehicle. The ideal absorber, from the standpoint of specific 
impulse, is hydrogen, but it is not absorbing to the visible and infrared radiation 
characteristic of current high power lasers unless it is heated to temperatures at 
which it becomes ionized. Since significant levels of ionization are not reached until 
the temperature has reached approximately 10,000 K, it is attractive to consider 
the use of mixtures of hydrogen and molecular absorbers, such as water vapor, to 
reduce the temperature at which the laser beam can be absorbed in the propellent. 

Both the required size of the absorbing chamber and the specific impulse of 
the propulsion device depend critically on the magnitude of the laser absorption 
coefficient for the mixture at the high temperatures desired for a rocket propellent. 
Prediction of the laser absorption coefficient is difficult for mixtures of water and 
hydrogen at high temperature due to dissociation of the water molecule under the 
influence of laser absorption and the complex chemical interactions involved. For 
this reason, several previous investigations [8,9,10] have attempted to measure the 
absorption coefficient of water at the 10.6 micrometer wavelength characteristic of 
the carbon dioxide laser. 

Two methods have been used to obtain the conditions of pressure and temper- 
ature required: heating of the mixture in a laser sustained plasma [9], and shock 
heating of the mixture in a shock tube [8,10]. 

Fowler [9] obtained measurements of the absorption coefficient at 10.6 microm- 
eters for mixtures containing H^O, D 2 O and NHq in hydrogen at pressures to 
13 atm. A laser sustained plasma was created in the test chamber using a vertical 
beam from a cw carbon dioxide laser of 7 kW power. Measurements of both tem- 
perature and absorption coefficient were made using Abel inversion of the output 
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of a Mach-Zehnder interferometer illuminated with a separate carbon dioxide laser. 
The values of absorption coefficient measured in this experiment (approximately 1. 
cm' 1 ) exceeded the calculated values by more than a factor of ten for both H 2 0 
and D 2 0. 

Kemp and Krech [8] made direct absorption measurements in the incident 
and reflected regions of a shock tube operating in, mixtures of H 2 0 and argon at 
temperatures to 3064 K and pressures to 49 atm. The maximum value of absorption 
coefficient measured was 0.45 cm -1 ama~ 1 . This is smaller, by a factor of 25, than 
the comparable value measured by Fowler . Rosen et. al. [10] made additional shock 
tube measurements in mixtures of C0 2 , H 2 0 and NHz with similar differences in 
measured values compared to those obtained by Fowler. 

To avoid the complications associated with the use of laser sustained plasmas, 
a method utilizing a laser absorption wave has been developed. In this technique 
a laser beam impinges on a ceramic surface, heating the surface and the cold pro- 
pellent mixture adjacent to it. When the mixture near the surface has reached a 
temperature where it becomes absorbing to the incident laser beam (approximately 
1000 K) , an absorption wave is initiated. This absorption wave propagates through 
the propellent mixture in the direction of the incident beam at moderate veloci- 
ties, and simultaneous measurements of temperature and absorption coefficient can 
be made as the wave passes through a test chamber. For this investigation, non- 
intrusive optical diagnostic techniques were developed to measure the translational 
temperature of the propellent mixture, and a test cell was designed and constructed 
to carry out the measurements. 


n. EXPERIMENT DESCRIPTION 


A. Absorption Wave Concept 

The experimental concept was arrived at jointly by personnel at NASA Mar- 
shall Space Flight Center (MSFC) and The University of Tennessee Space Institute 
(UTSI). It consists of using a laser absorption wave initiated from a plate heated 
by the impinging laser beam to create the conditions of temperature and pres- 
sure required for the absorption measurements. As the absorption wave propagates 
through a mixture of water vapor and a diluent gas, the mixture is heated to the 
maximum temperature which can be sustained by the laser absorption. When the 
wave propagates past observation ports in the test chamber, simultaneous measure- 
ments axe obtained for the optical transmission of a separate 10.6 micrometer laser 



and for the translational temperature of the mixture. Thus, for each occurrence 
of a wave, a series of transmission measurements are obtained at different values 
of mixture temperature. A schematic of the experimental apparatus is shown in 
Figure 1. 


Since the mixture will be heated by strong absorption of the laser beam in the 
rovibrational bands of the H^O molecule, the temperature within the absorbing 
beam will not be uniform, and it is likely that the mixture will not be in thermody- 
namic equilibrium. To obtain accurate measurements of absorption it is necessary 
that the transmission measurements be made across the entire cross-section of the 
wave. A pyroelectric array detector was chosen to obtain the required spatially 
resolved measurements of laser transmission. The transmission data can then be 
Abel inverted to give spatially resolved measurements of the absorption coefficient. 
An optical measurement of temperature utilizing the measurement of acoustic wave 
velocity was chosen to avoid the uncertainties in rotational or vibrational tempera- 
ture measurements arising from the uncertain state of equilibrium in the mixture. 
Acoustic wave speed is determined from a time-of-flight measurement of the decay- 
ing blast wave produced by an excimer laser induced spark, using helium-neon laser 
beams of known separation. 

A simple one-dimensional model of the initiation and propagation of a laser 
absorption wave was developed to determine whether this concept was feasible. 
The primary goals of this analysis were to determine if the wave speed was slow' 
enough that the required measurements could be obtained during wave passage, 
and to determine the threshold laser beam intensity for the initiation of the laser 
absorption wave. 


B. Numerical Model 

A numerical model was constructed in order to study the feasibility of the 
laser absorption wave technique and to provide a first approximation for the time 
and spatial characteristics of the phenomenon. The model was based on a one- 
dimensional, time-dependent algorithm for thermal conduction in the target gas 
with volumetric absorption of laser beam power governed by Beer’s Law. 

The derivation of the finite element approximation follows the method given 
by Patankar [11]. Consider the one-dimensional unsteady heat conduction equation 
with source term: 
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( 1 ) 



where T is temperature, t is time, x is the spatial coordinate, p is the density, c is 
the specific heat and k is the thermal conductivity. 

For the case of an absorbing gas with emission of radiation, the source term is 
given by: 


S = al - $ (2) 


where a. is the absorption coefficient for laser radiation of intensity I and $ is the 
net radiation emission coefficient. 

The intensity of incident laser radiation is prescribed by Beer’s Law: 


( 3 ) 


For a gaseous medium the properties are related by the equation of state: 


"-«r 

dx 


for a laser beam propagating in the negative x direction. 


P = pRT 


( 4 ) 


where P is pressure and R is the gas constant. 


The heat conduction equation may be discretized by integration over a finite 
time interval and over the control volume of Figure 2, 
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rt+At re dT rt+At per Q / dT 

L L pc m ixit = J, y. 


H 


dxdt. 


( 5 ) 


In the equations below, denote the state at time t by the superscript 0, and at 
time, t + At by the superscript 1. Assume the quantity pc is constant over the time 
interval At. Then equation (5) may be written as: 


- ir l ( K f ) + s \ ixit - 


( 6 ) 


The integral over x in (6) may be evaluated and separated into terms to yield: 


/ t+At ■ 

K. 


(Tk - T') (T* - r*) 


Sx. 


- /c. 


Sx„ 


+ SAX 


dt (7) 


Here, the partial derivatives have been replaced by the difference approxima- 
tion, and the source term has been replaced by the mean value over the cell, S. If 
it is assumed that all properties except the grid temperature T p are constant over 
the time interval At then Eq. (7) may be integrated over time to obtain the fully 
implicit difference equation: 


(H-KJ-SA.* + 


pcAx (T 1 - T 0 ' - 5 
At V p p 


6x, 


Sx 


( 8 ) 


XV 


Approximate the source term S by the sum of a constant and linear term as 
follows: 


S — S c + SpTp 


( 9 ) 
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Equation (8) may now be rewritten by the collection of terms into a new form 
as follows by the definition of a series of grid coefficients: 


ApTp = AeTe + A W T W + b 


( 10 ) 


where 


Ae = K e /Sx e 


A\y — k w j 6x w 


0 _ pcAx 

Ap ~ ~aT 

Ap = Ap + Ae + Aw — S p Ax 


and 

6 = A° p T° + ScAx. 


This new formulation simplifies the numerical solution of the problem as it 
allows the matrix formulation of the equation set to be cast into a tri-diagonal 
matrix form. Note that the behavior of temperature at the given point depends 
only on the values to the left and right in the spatial grid. This formulation will 
create a ” banded” matrix of coefficients. The solution algorithm proceeds at each 
time step by the solution of a spatially coupled, steady-state problem using a tri- 
diagonal matrix algorithm. The parameters and properties used in the formulation 
of each steady-state problem axe the results of the solution from the previous time 
step. In this manner, the temporal dependence is determined. Note, however, that 
the properties ”lag behind” the solution by one time step. This requires the use of 
very small time steps for the solution to accurately track the true time behavior. 
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The source term has also been linearized by the assumption that the tempera- 
ture dependent part S is zero. Instead, the constant portion of the source term, S 
is calculated from the properties belonging to the previous time step. 

The spatial grid was divided into two principle regions: one of an absorbing gas 
and the second of a ceramic target as shown in Figure 3. The boundary conditions 
were chosen so that the temperatures on either end of the grid were prescribed for 
all times t. The initial conditions were those of uniform temperature equal to the 
temperature prescribed on either end for all points in both regions. The laser beam 
is assumed to proceed from right to left in the solution grid as shown in Figure 3. 
The source term S may be written in discrete form by an approximation to Eq. 
(2) . The energy absorbed into the nth volume element is given by: 


__ hf+i — In 

Axn 


( 11 ) 


The laser intensity I is determined by integration of the Beer’s law equation at the 
end of each time step. The absorption coefficient for 10.6 micrometer radiation used 
in the model was curve fitted to shock tube data from Kemp and Krech [8] for the 
P(20) CO 2 laser line. 


The net source term is the difference of the energy absorbed and the energy 
lost due to radiation $. The total radiation from the nth volume element is given 
by: 


is = 4 < K„ > 


( 12 ) 


Here < K v > is the spectrally averaged absorption coefficient . The absorption 
coefficient is related to the optical cross section <c v > by: 


< K v >= N a <o v > (13) 

where N a is the number density of the absorbing species. 

Radiation from both H^O and OH were considered in the program in order to 
calculate the net source term. The cross-sections were calculated from experimental 
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data obtained by Kemp and Krech [8] who give values for the emissivity of these 
species for various values of pressure-length product, ps. The cross section < a v > 
is then related to the emissivities by: 


< a v >= 


eKbT 

ps 


(14) 


The thermal conductivities and specific heats for the gaseous H 2 /H 2 O mixture 
were taken to be equal to those of the major constituent if 2 * These data were taken 
from Patch [12] who calculated properties for spin-equilibrated H 2 for a range of 
pressures and temperatures. 

The properties were treated differently in the ceramic region of the solution 
grid. This was accomplished by the use of steering logic in the property subroutines 
. The solution algorithm is not disturbed by discontinuity in properties along the 
solution grid as long as the grid is sufficiently fine in the region of the discontinuity. 
The "fineness” of the grid at the discontinuity was chosen by trial and error to 
achieve stability in the solution without incurring excessive memory requirements. 


The source terms for this region of the solution were set to zero except for the 
first ceramic cell at the solid-gas interface. This cell absorbed all of the remain- 
ing laser power and emitted radiation given by the Stefan-Boltzmann relation for 
solid body radiation from a surface. The emissivity of the ceramic was chosen to 
correspond to 100% dense polycrystalline AI 2 O 3 at elevated temperatures, as were 
the specific heats and thermal conductivities in this region of the grid. None of the 
re-radiated power from the ceramic surface was considered to be re-absorbed in the 
adjoining gas mixture. 


The mole-fraction of each constituent of the gaseous mixture will change as 
the temperature is elevated. This is due to the dissociation of the water vapor 
and diatomic hydrogen to produce other species such as 0 H, 02 ,H,H + , etc... . 
This effect was included in the program by curve fitting the output of the NASA 
code ODE which was rim for each case to be considered. It is the inclusion of 
these effects which limits the maximum temperature of the thermal wave since 
the number density of absorbing molecules is reduced by the onset of dissociation 
at higher temperatures. In this model, only H 2 0 is considered as an absorber, 
hence the maximum temperature is limited to approximately 4000 K for the cases 
considered. 
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Two major cases were considered for use of the model. These were for a mixture 
of Hz and HzO, with 20 percent H 2 0 by volume initially at ambient temperature 
for pressures of 1 and 10 atm. These cases represented extremes readily accessible 
to laboratory apparatus and encompassed regimes of interest to propulsion. Two 
major sets of property subroutines were built for these cases, primarily by simple 
curve fits to available data. Conservative estimates were used where data were 
sketchy or unavailable. 

The results of the two major cases are shown plotted in Figures 4 and 5. The 
assumed incident laser intensity was 375 W/cm 2 for both of these cases. The 
time steps were .01 sec for the 10 atm case and .05 sec for the 1 atm case. The 
temperature profiles are plotted every 2 sec up to a maximum of 20 sec, in order 
to show the progression of the wave with time. For the 1 atm case shown in Figure 
4, the wave velocity reaches a constant velocity of 7.5 cm/sec after 10 sec. The 
maximum temperature in the wave is about 3600 K and the maximum temperature 
on the ceramic wadi reaches 1800 K. For the 10 atm case, the wave velocity reaches 
3 cm/ sec after 10 sec. The maximum temperature in the thermal wave approaches 
4000 K and the maximum wall temperature is about 1600 K. The higher wave 
velocity for the 1 atm case is a result of the lower thermal conductivity which leads 
to steeper thermal gradients at the leading edge. The higher wall temperature for 
the 1 atm case is also due to the lower thermal conductivity which reduces heat 
transfer from the wall surface. 

The laser intensity has a major effect on the behavior of the solution. Several 
alternate cases were examined to investigate the solution for lower intensities. One 
of these cases is shown plotted in Figure 5. No persistent wave is launched from the 
wall for this case. Instead, a heated zone is established away from the wall which 
converges to a steady-state profile. This behavior is due to a lack of incident power 
necessary to overcome the radiative and conduction losses. 

The conclusions drawn from this study were that it was possible to create a 
laser absorption wave and that it might provide a viable means for the creation of 
the conditions necessary for the study of H 2 /H 2 0 mixture properties at high tem- 
peratures. The study did demonstrate a sensitivity to laser power levels, although 
the predicted power levels were within the reach of available high energy lasers. 

It is important to note that several important loss mechanisms were not in- 
cluded in the model. The major factor not included in the study was that of the 
convective term in the energy equation. The exact effect of the convective terms 
have not been determined due to the difficult nature of such a solution. Unfortu- 
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nately, time did not permit a more complete analysis of the convective flow through 
the wave. It was anticipated that this effect would raise the laser power threshold 
for the initiation of the wave due to the losses incurred and the subsequent reduc- 
tion in thermal gradients. Another important effect not included in the model was 
that of radial conduction in the thermal wave and the ceramic end-wall. This effect 
would also raise the necessary laser power threshold but it was felt that this effect 
would be minimal because of the size of the beam and the nature of its unstable 
resonator beam pattern. 

Complete Fortran 77 program listings for the computer program axe included 
in Appendix A. This program was executed on a VAX 11/780 machine at the 
University of Tennessee Space Institute. Property subroutines are also included in 
the listings for the 1 and 10 atm cases. 

C. Temperature Measurement 

The measurement of local, time-resolved gas temperatures was essential for 
the prediction of the varying constituency of the high-temperature, subsonic, gas 
sample. Since the laser absorption measurements were to be performed in the com- 
paratively slowly varying trailing-edge of the thermal wave, the gas mixture param- 
eters and their gradients in this trailing-edge region axe of fundamental importance 
in deter min ing the criteria that must be satisfied by the temperature measure- 
ment method. Figure 7 shows a sample result for the space-time variation of the 
thermal wave for a 20 %H 20 / 80 %Se volume mixture at 10 atm total pressure for 
I 0 = 375 w/cm 2 . At a sample axial location, say x = 50 cm, from Figure 7 it is seen 
that an axial spatial resolution of approximately 5-10 mm will yield an uncertainty 
of less than pm ± 5% in T, which is an acceptable measurement error. Further, 
from Figure 7 it can be estimated that the axial flow speed is less than 10 cm/ see. 
Consequently, the temporal resolution necessary to determine accurately the ther- 
mal wave profile in the trailing-edge region is seen to be on the order of 10 msec 
or less. An additional criterion that should be satisfied for measurement of the 
temperature of the reacting gas sample is an insensitivity to the gas sample’s local, 
instantaneous constituency, which should be considered as an a priori unknown for 
the experimental measurements. 

A variety of methods of temperature measurement exists; most of these meth- 
ods, however, fail to satisfy the desired criteria. For example, since local thermal 
equilibrium cannot be ensured and the spectral broadening parameters are poorly 
known for the high-temperature region of this study, emission-absorption techniques 
cannot be used. Spontaneous Raman scattering (SpRS) is attractive, but the SpRS 
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cross section and scattering signal are so low that adequate time resoltuion can- 
not be achieved for a single thermal wave; i.e., multiple thermal waves would be 
required to map the spatial profile of an "average thermal wave” profile. Coherent 
Anti-Stokes Raman scattering (CARS) is known to be capable of temperature mea- 
surements with good signal-to-noise (S/N) ratios for combustion sources. However, 
for this study, the pressure-temperature region is such that significant thermal dis- 
sociation occurs, resulting in a substantial decrease in all initial molecular species. 
The OH radical species will, however, increase during the disociation processes, 
but its concentration is too small for its use for a CARS temperature measurement 
with large values of S/N. Furthermore, LTE may be a questionable assumption for 
OH radicals produced by thermal dissociation of H 2 0 that has been vib rationally 
excited by absorption of C0 2 laser radiation. As a result, CARS was rejected as 
a practical approach for the temperature measurement. Laser-induced fluorescence 
(LIF) satisfies all criteria but the last; the results may be sensitive to the instan- 
taneous mixture constituency. Before embarking on a study related to quenching 
cross sections of various molecular species and the effects of such cross sections on 
temperature measurements, other possible methods were sought. 

The method of Zapka et al [14] was chosen for the temperature measurement 
of and is designated the optoacoustic laser-beam deflection (OLD) method. As 
sketched in Figure 8, OLD consists of using a high-power laser which is focused to 
a small volume in the gas sample. The excitation laser irradiance exceeds the gas 
breakdown threshold, and a plasma spark and cylindrical blast-wave are produced. 
The blast-wave propagates outward and, after a distance ip, it decays into an 
acoustic wave. Two probe-laser beams, separated by a distance Ax, detect the 
passage of the acoustic wave with a passage time difference At, which corresponds 
to the separation distance Ax . If the gas sample possesses a flow speed vp and a 0 
is the sample’s sound speed, then 


a 0 -\-vp = Axd/Atd 


(IS) 


and 


Qf0 Vp Aty 


(16) 


where the subscripts d and u denote the downstream and upstream directions, 
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respectively. 

For an assessment of the order of magnitudes (O), 


O(a 0 ) >~ 10 4 cm /sec. 


0(v p) «~ 10 2 cm/ sec , 

so that 

0(vf/ci>o) <~ 10 — 2 , 

and now 

do — Ax<i,u/ (17) 

For the subsonic flow of assumed one-dimensional speed vp, note that Eqs. (15 
and 16) indicate a method of measurement of Vp: 


vp — (l/2) [(Axoi/ Atc{) (Ax u /At u )]. (IS) 


Consequently, probe beam-pair detectors located both downstream and upstream 
of the excitation beam can be used to infer values of vp. 

Now, returning to Eq. (17), it is recalled that for a perfect gas of specific heat 
ratio 7 , molecular weight M and temperature T, 


a 2 = ( dp/dp)S = iRT/M (19) 
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where p and p denote, respectively, the gas pressure and density, and subscript S 
denotes an isentropic derivative and R is the universal gas constant. Therefore, if 7 
and M of the ideal gas are known, measurement of a 0 yields a measurement of T. 

Several features of the gas sample require much closer scrutiny before declaring 
as feasible the use of OLD for measurement of temperature in this study. First, the 
gas sample is reacting, and an initial assessment of OLD would indicate the likeli- 
hood of dependence of a 0 on the sample’s constituency. Second, the use of a pulsed 
excitation source of short-duration (= nsec ) will produce a short-duration acoustic 
pulse, and the Fourier transform of this pulse will exhibit a range of frequencies. 
Since the sample consists of molecular species and acoustical dispersion will exist, 
the validity of describing the sound speed by an isentropic (and dispersion-free) re- 
lation (Eq. (19)) must be determined. Third, the question of the effect of refraction 
on the excitation and probe beams must be considered, and, fourth, the magnitude 
of A xd and the limiting spatial resolution that is set by (Ax<f )U ) moz and (A ti lU ) mtn 
axe of interest. 

To determine the sensitivity of a 0 on constituency, it was necessary to calculate 
the equilibrium concentrations of the dissociating, reactive samples as a function 
of temperature. For this purpose and subsequent rate kinetics calculations, the 
reaction rate data of Refs. [15] and [16] were used. Since H 2 0/He, H 2 0/Ar and 
H 2 O/H 2 mixtures were all of interest, equilibrium concentration and sound speed 
calculations were performed for all three mixtures for a range of initial concen- 
trations and pressures for 1000 ~< T <~ 5000 K. Figures 9 through 12 show 
sample concentration results for Ar and H 2 mixtures. Specifically, for the ini- 
tial mixture (by volume) of 20%H 2 O and 80% Ar at P = 1 atm. Figure 9 shows 
the variation of the concentrations of H 2 0, H 2 ,. H, and OH, [H 2 0], [H 2 ], [ H ] 
and [OH], respectively, as a function of T. Figure 10 shows similar variations for 
P = 10 atm. Several features of Figure 9 and 10 are of interest. First, note the 
small concentrations of H 2 0 for both P — 1 and 10 atm at T = 4000 K. Second, at 
T = 4000 K, [H 2 ] <~ .03 and note the increasing concentration of H([H] >~, 0 . 1 ) 
for T = 4000 K. Finally, the peak in [H] at a value of approximately 0.04 is of 
interest and justifies the rejection of CARS as a candidate measurement method. 
Further, a comparison of Figures 9 and 10 shows clearly the suppression of the 
dissociation reactions that results from the incresae of pressure from 1 to 10 atm. 
Figures 11 and 12 show similar variations for a 20%fT 2 0/80%ff 2 initial concentra- 
tion mixture at P =1 and 10 atm, respectively. Relative to the Ar mixtures, a 
significant increase in [H 2 O], [H 2 ], [OH] and [FT] understandably results. Further, 
the peak in [OH] is shifted to T/ 4000 K for P — 1 atm and significantly beyond 
4000 K for P = 10 atm. Finally, from such results, the equilibrium sound speed ao 
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was computed, and, as examples, Figures 13 and 14 show the variation of ao with 
T for 20 %H 2 O /80%H e at P = 5 and 10 atm, and Figure 15 shows ao versus T 
for a 20%H2O/80%H2 mixture at P = 1 atm. From Figures 13 and 14, it can be 
seen that a measurement of the equilibrium sound speed provides a unique value of 
temperature for an equilibrium reacting mixture. Further, note that little difference 
in a 0 versus T results from the change of pressure from 5 to 10 atm for H 2 O/ He 
mixtures. As expected, however, significant variations result for ao versus T as one 
changes the non-absorbing gas from He to H 2 . Results similar to Figures 13 and 
14 were found for H20/Ar mixtures. 

The calculation of the equilibrium, isentropic sound speeds relied upon the va- 
lidity of Eq. (18) which presumes the absence of dispersion effects. Specifically, 
Eq. (18) assumes the absence of relaxation processes and be categorized as a zero- 
frequency speed of sound. As the sound wave frequency (u) increases, relaxation 
processes result and the internal modes of molecular motion fail to maintain equilib- 
rium with the translational mode. For an example of order of magnitude relations, 
suppose the molecular species is characterized by an equilibrium value of 7 = 1 . 2 . 
As v increases, 7 increases from 1.2 to 1.667, the monatomic gas value, and the speed 
of sound increases according to Eq. (18). The inferred values of corresponding to 
this range of 7 and a% differ by almost forty percent. Therefore, it is necessary to 
determine the Fourier frequency spectra expected for a laser excitation pulse-width 
on the order of 1 — 10 nsec. 

Assume the laser p- induced, time-dependent pulse amplitude can. be described 
by the functional form 


El = Eocos^t/ti). 


( 20 ) 


The Fourier transform El (w) is easily found to be 


El{u) = 2 t l E 0 - cos(wt L /2)/[l — w 2 ti/fl- 2 )] (21) 

where u> = 2nv and Tl = 2 ti/ 7 r. 

From Eq. ( 21 ), it is seen that the predominant Fourier amplitude occurs at 
u = w = 0 and that the first and largest "secondary maximum”, which occurs 
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at U\tjj = 57 r/ 8 , has an amplitude less than 10 -1 times the zero frequency am- 
plitude. Since fx, ~ 10 nsec = 10 -8 sec, Vi = 200 MHz. For P = 10 atm 
and T >~ 3000 K, at V\ = 200 MHz, all vibrational modes can be considered 
frozen and rotational relaxation will have become significant. [15]. As a result, at 
Ui = 200 M Hz, q ( 200 Af Hz) >~ 1.4 and vibrational dispersion would result in an 
increase in a \ (and inferred (and incorrect)) T of approximately twenty percent. 
Therefore, in summary, for a temporally symmetric 10 nsec laser excitation pulse 
and blast wave, the acoustic pulse can be expected to exhibit detectable dispersion. 
However, the predominant (and slowest) component of the anticipated asymmetric 
acoustic wave will be the zero-frequency component that is predicted by Eq. (18) 
and shown in the previous figures. 

The preceding conclusions were formed on the basis of the assumption of chem- 
ical equilibrium. Since, as Figure 7 shows, a laser absorption wave processes the ini- 
tial mixture sample, time-dependent reactions result and the chemical constituency 
will vary with time. The evaluation of the effect of such finite-rate chemical kinetics 
on the previous conclusion required an assessment of the time required to achieve a 
condition approximately that of local thermal equilibrium. To do this, differential 
equations corresponding to the reaction processes of Refs. [16] and [17] were solved 
numerically as a function of time to yield the time-dependent concentrations of 
the mixture. To describe approximately the passage of the laser absorption wave, 
the gas temperature of the initial mixture was assumed to be subjected to a step- 
function increase which corresponded to the arrival of the leading-edge of the wave. 
Recall that from Figure 7 and the estimated value of vjr(~ 10 cm/ sec), one expects 
characteristic times of significant change of parameters (r p ) in the thermal wave 
pulse to be no less than 10 msec. Therefore, if one finds the (maximum) time to 
approach chemical equilibrium (r c ([x]) max ) for the slowest-varying major specie x 
to satisfy the inequality 


r c ([x]) max ^ Tp, 


the assumption of LTE and the previous conclusions are valid for measurements 
made with data acquisition and resolution times tm on the order of, or slightly less 
than r p . 

Figures 16 through 34 show sample results for the variation with time of the 
mole fractions of the mixture species for H^O jH t and H 2 O / mixtures. 

Figures 16 and 17 show the 20% H 2 O / S0%H e results for a step-function increase 
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in T from 300 K to 2000 K for P = 10 atm. It is seen that all species except H 2 0 
and 0 2 attain a steady-state, and LTE, in times on the order of 100 /zsec; 1 msec, 
however, is required for H 2 0 and 0 2 . Even so, t c {[H 2 0] or [0 2 ]) max = 1 msec <C 
10 msec, and LTE cam be assumed and the desired temporal resolution achieved if 


1 msec < tm < r p — 10 msec. 


Figures 18 and 19 show the results at P = 5 atm and T = 3000 K, and Figures 20 
and 21 present the P — 10 atm and T = 3000 K case. For P = 5 atm, T c ([x]) max = 
1 msec but for P = 10 atm, where the dissociation reactions are retarded, or slowed, 
T c{[ x ])max > 1 msec. Therefore, for P = 10 atm and T = 3000 K, more and careful 
analysis would have been required for accurate analysis of the data. Similar results 
are shown for T = 4000 K in Figures 22 through 26 as P increases from 1 to 
10 atm. The increase in the rates of the reaction kinetics at 4000 K results in 
T c{[x])maz < 1 msec. Therefore, 


1 msec < tm < r p = 10 msec. 


The results for mixtures of H 2 0 and H 2 are shown in Figures 27 through 34 
and inequalities similar to the immediately preceding value result. 

Finally, for the condition represented by Figure 34, the time-dependent sound 
speed was computed throughout the reaction processes. Figure 35 shows the sound 
speed versus time for a step-function temperature change. From Figure 35, it is to 
be noted that, following the step-function increase in T and initial, rapid change 
in constituency, neglect of the time-dependent chemical kinetics could result in an 
error in the sound speed of about 8%, which translates into an error of 16% in 
temperature. Consequently, for the range of conditions that were anticipated for 
this study, it is essential to perform such time-dependent kinetics calculations to 
either avoid or correct for finite-rate chemical kinetics effects on the sound speed 
and temperature measurements. 

On the basis of the preceding studies and calculations, the OLD technique 
was adjudged to be feasible for determination of temperature for the dissociating, 
reacting sample of H 2 0/He, H 2 0/Ar and H 2 0/H 2 mixtures at the composition 
and pressure conditions anticipated for use. 

The final major concern regarding application of the OLD method, as well 
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as other optical techniques, is the effect of the gas dynamic parameters on the 
quality of the pump and especially the probe beams. Of particular concern, of 
course, is the effect of refraction and resulting deviation of the pump and probe 
beams that result from the gradient of the index of refraction of a laminar flow- 
field as well as the refraction and beam spreading encountered in a turbulent flow- 
field. For the conditions anticipated for this study, it is seen that O(Rcl) = 10 2 , 
where Re^ is the unit length Reynolds number of the gas. Therefore, since a 
characteristic length L is on the order of 10 cm, it is quite possible that the flow-field 
is laminar. If so, for a cylindrically symmetric CO 2 laser beam injected along the 
centerline of the gas sample, a detailed analysis of the experimental results requires 
the use of the eikonal equation for determination of the refractive displacement of 
the beams. The approach used by Lewis and Selman [18] provides a prediction 
of the refractive displacement of the irradiance centroid of the probe beams for 
specified refractive index (or gas density) gradients. From the expected geometry 
of the experi mental configuration, it is readily seen that diametral injection of the 
probe beams of minimal or very small diameters minimizes beam refraction effects. 
Further, quantitative evaluation of such refraction effects is easily accomplished by 
observing probe beam deflections without the use of the excitation pump beam. 

D. Absorption Measurement 

Measurement of the absorption coefficient at a particular point in a volume 
having axial symmetry can be accomplished if a transverse scan of transmission can 
be obtained. If the measured transmission transverse to the axis is given by T(x), 
then the absorption at the radial distance, r, from the axis is given by the Abel 
inversion, 



<n[r(x)l 
(x 2 — r 2 ) 1 / 2 


dx 


(23) 


where ro is the distance from the axis where the absorption coefficient k(r) vanishes. 
Evaluation of this expression for experimentally measured transmission scans can 
be performed using a technique developed by Shelby and Limbaugh [19]. 

The transverse scan of transmission was obtained using a linear pyroelectric 
array having 32 separate elements. A stabilized, continuous, 6 W carbon dioxide 
laser beam passed through a variable attenuator, then it was expanded and colli- 
mated before passing across the test chamber through zinc selenide windows. After 
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it exited the chamber, the beam passed through a rotating chopper wheel and was 
imaged onto the detector array so that the elements of the array spanned the region 
of interest. The output beam from the laser could be directed to a small grating 
monochrometer to determine the particular rotational transition which was lasing. 

The output of the detector array was sampled at a frequency of 5 Hz and read 
into an LSI-11 data acquisition system which was also used to acquire the data 
from the temperature measurement. Ti ming for the experiment was controlled by 
the computer using the trigger signal from the rotating chopper, to insure that the 
temperature and transmission measurements coincided in time. 

E. Test Chamber 

The test chamber and the associated boiler designed to supply the chamber 
with a known mixture of water vapor and diluent gas was designed and fabricated 
at MSFC and provided to UTSI for the preliminary experiments performed there. 
The chamber was constructed from two existing sections of flanged aluminum pipe 
which were attached to a machined aluminum hexagonal section fitted with six 
windows to allow optical access. Four of the windows were quartz to permit access 
for the excimer laser and for the helium-neon laser beams used for the temperature 
measurement. The other two windows were made of zinc selenide to permit access 
for the 10.6 micrometer carbon dioxide laser used for the transmission measurement. 
A third zinc selenide window was fitted to the flanged cylindrical section to permit 
access for the high power 10.6 micrometer laser beam used to create the laser 
absorption wave. The test chamber is shown schematically in Figure 1. 

A boiler system was designed to supply the test chamber with a known mixture 
of water vapor and diluent gas. This system was composed of a spherical stainless 
steel tank of approximately 10 gallons internal volume wrapped with heating coils 
and heavily insulated. Fittings were installed in the top and bottom of the tank 
to allow for injection of the diluent gas. The tank was operated approximately 
half full of water and the diluent gas was introduced through a 1/4 inch stainless 
steel standpipe and bubbled through the water in the tank. The resulting water- 
saturated vapor was extracted from the top of the tank through a glass-wool filter 
to eliminate any water droplets. The pressure in the entire system was controlled 
by throttling the flow at the outlet of the test chamber. 

The temperature of the boiler was controlled to an accuracy of 0.5C using a 
platinum RTD thermocouple and controller. The partial pressure of water vapor 
at the boiler operating temperature was determined from the equilibrium steam 
tables. The pressure of the diluent gas on the water was increased to provide the 
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desired partial pressure ratio (mole fraction) of diluent gas to water vapor. Using 
this procedure, the mole fraction of diluent gas exiting the boiler was always known 
during the test. The boiler typically was operated at a temperature of 120 C and 
a pressure of 75 psig. This condition corresponds to a water vapor mole fraction of 
0.33. 


The test chamber was connected to the boiler by a 6 ft length of 1/2 inch di- 
ameter stainless steel tubing. The pressure drop through this length of connecting 
tubing at the relatively low flow rates (approximately 2 liters/min) was negligible. 
All tubing, fittings and test chamber surfaces were heated to a temperature 30 C 
higher than the boiler temperature in order to eliminate any possibility of vapor 
loss due to condensation. Initial testing of the system revealed that thermal strati- 
fication was a problem in the boiler. This problem was eliminated by the addition 
of a closed-loop pumping system which removed water from the bottom of the tank 
(past the RTD probe) and replaced it at the top. The pumping rate was 4 gpm 
so that the volume of the water was circulated approximately once every minute. 
The pumping loop was insulated in order to avoid cooling of the water during cir- 
culation. The state of equilibrium in the tank was assured by the installation of a 
bare wire thermocouple in the top of the tank to measure the temperature of the 
gas/water vapor mixture. When the gas temperature and the water temperature 
agreed to within the accuracy of the temperature probes, the state of equilibrium 
was assumed to exist. This state of equilibrium was normally reached within a 
period of 4 hours from initial startup of the boiler/test chamber system. 

F. Data Acquisition and System Integration 

Successful measurement of the absorption coefficient as a function of tempera- 
ture relies on the simultaneous acquisition of the transmission data from the probe 
laser and the time-of-flight measurement for the acoustical wave. This requires syn- 
chronization between the firing of the excimer laser and the chopper wheel used 
for the transmission measurement, using the data acquisition computer to control 
the timing for the acquisition of the pyroelectric array scans and the initiation of 
the counters used for the time-of-flight measurement. The trigger pulse from the 
rotating chopper wheel was used to initiate the timing sequence within the LSI-11 
computer which delivered the firing pulse to the excimer laser, started the start and 
stop counters for the time-of-flight measurement and acquired the scan from the 
pyroelectric array. When one measurement sequence was complete, the counters 
were reset and a new sequence initiated at the next trigger pulse from the chopper. 

A measurement begins when the high power laser is turned on and begins 
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heating the ceramic end-wall of the chamber. All of the scans from the pyroelectric 
array are recorded by the computer, but prior to the time that the wave is launched 
from the wall the absorption across the chamber is negligible, and these early scans 
are averaged and used to normalize the later scans to obtain the transmission. 
Since the temperature is measured at only a single point within the wave, the 
transmission must be Abel inverted using the method described in section II-D to 
relate the absorption coefficient to the measured temperature. 


HI. EXPERIMENTAL RESULTS 


An experimental study was performed at UTSI to verify the proposed method 
of temperature measurement, to evaluate the utility of some candidate materials 
for the end wall and to attempt to launch a laser absorption wave using the small 
(1.5 kW ) laser at UTSI. 

An experimental study of the temperature measurement technique was initiated 
to determine if the excimer laser was an acceptable source for the breakdown spark, 
and to determine the suitability of the available instrumentation for measurement 
of the time-of-flight of the resulting acoustic wave. For the initial investigation, 
the experiment was set up in ambient air using a boxcar integrator to measure the 
passage of the wave through two helium-neon laser beams which were separated by 
approximately 2 mm. It was found that the excimer was effective at producing the 
laser sparks with repetition rates to 200 Hz, and that the measured temperature 
agreed with the known temperature of the ambient air to within 1%. 

Following the initial experiments in air, a small (10 cm diameter, 30 cm long) 
chamber was constructed to contain an atmosphere of argon or helium, two of the 
proposed diluent gases for the absorption measurement. For these experiments the 
time-of-flight of the wave was measured both by the boxcar integrator and by using 
time-to-digital counters integrated into the LSI-11 computer. In order to insure 
reliable triggering of the counters by the deflection of the helium-neon beams, the 
beams had to be accurately focused on the detector pinholes. The chamber could 
be heated by use of a portable, hand-held, heat gun and the chamber temperature 
measured with a thermocouple. In argon, the experiments verified that tempera- 
tures could be determined to approximately 1%, and were independent of chamber 
pressure over the range from one to five atmospheres. However, when helium was 
used as the test gas, the excimer laser failed to produce a spark breakdown, even 
at helium pressures up to five atmospheres. Addition of water to the cell was found 
ineffective in promoting breakdown in helium. 
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An ideal material for the end wall would have the following characteristics: it 
would absorb all of the incident laser energy, have low thermal conduction to reduce 
the energy conducted away from the surface, withstand thermal shock without 
fracturing, and not melt at the temperatures required for initiation of the absorption 
wave. Experiments were conducted at UTSI using the 1.5 kW carbon dioxide laser 
to screen some possible candidate materials. 

Alumina ( AI 2 O 3 ) has a high absorptivity at 10.6 micrometers, but thin sections 
were found to shatter when subjected to the laser beam, and thicker sections were 
found to melt and drip. Several varieties were tested and none were found suitable. 
Tungsten had too great a reflectivity at 10.6 micrometers to be effectively heated 
by the incident power available from the 1.5 kW laser, and it did not get sufficiently 
hot to obtain an optical pyrometer measurement of surface temperature. Graphite 
absorbs the laser beam efficiently, but the energy was carried away rapidly by 
thermal conduction and, again, it was not possible to obtain a measurement of 
surface temperature using the optical pyrometer. A ceramic material of unknown 
composition, Norton Magnorite was found to have the properties required. In a 
section of 1/2 inch thickness it was possible to obtain surface temperatures of 1850 
C, as measured with the optical pyrometer, and it was possible to heat it for as long 
as six minutes without melting or thermal fracture. After repeated heating cycles 
the surface appearance changed from the initial brown color to a crystalline white, 
and some cracking was noted after cooling. Since the measured surface temperature 
of this material was sufficiently high, according to the model predictions, to initiate 
a laser absorption wave, it was chosen as the end wall material for the absorption 
wave experiments. 

The model predictions indicated that approximately 375 W /cm 2 would be re- 
quired to initiate a laser absorption wave, and the 1.5 kW laser at UTSI was unable 
to produce this flux over a spot size sufficient to obtain a reliable measurement of 
transmission. However, the high power laser at MSFC was under repair and not 
available, and it was decided to attempt to initiate a laser absorption wave using the 
UTSI laser. This would provide an integrated test of the data acquisition system 
and the test chamber, even if wave initiation was not achieved. 

The beam from the 1.5 kW laser was focused by a spherical mirror having 
a focal length of 10 m to a spot approximately 2 cm diameter. Burn patterns 
in plastic showed that the beam was fairly uniform in intensity across the spot, 
due to the effects of spherical aberration and coma, but showed some striations 
due to diffraction. The power of the beam after passing through the optical train, 
including a wedge beamsplitter used to monitor the power during the experiment, 
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was measured by a Scientech calorimeter and found to be 750 W. Thus, the incident 
laser intensity on the ceramic end-wall was approximately 240 W /cm?, which was 
marginal according to the model predictions. For this spot size, using the maximum 
absorption coefficient reported in Reference [8], the minimum transmission would 
always exceed 92%, but that value of transmission was thought to be detectable. 

When the experiment was first attempted, it was found that the laser beam used 
for the transmission measurement was not uniform after passing through the test 
cell. This proved to be due to the effect of interference of the coherent laser beam 
passing through the nearly-plane zinc selenide windows. To alleviate this problem, 
new windows having a wedge angle of 1 degree were obtained. The interference 
fringes from the wedged windows were sufficiently close together that approximately 
three fringes fell on each of the detector elements of the array, effectively smoothing 
the beam profile. 

Two attempts were made to initiate the laser absorption wave: one in a' mixture 
of 0.33 mole fraction of water in argon at a total pressure of 6 atm , and one in a 
mixture of 0.33 mole fraction of water in helium at a total pressure of 6 atm.' 
During these experiments, scans from the array detector were obtained at a rate 
of 5 Hz for approximately 1.5 min. Prior to a test, 1000 scans of the array was 
recorded to provide a background profile of the laser beam used in the transmission 
measurement. 

Figure 37 shows an average of these 1000 scans where the abcissa is analog-to- 
digital converter counts. The laser beam intensity was set to provide a maximum 
signal level of approximately 90% of full scale. The beam profile is fairly uniform 
in the region of the incident laser beam where the absorption is expected. The 
probe beam was focused on the array with a magnification chosen so the full array 
covered a distance in the test cell of 4 cm. Figures 38 and 39 show comparable ar- 
ray scans taken during the argon and helium tests, respectively. The scans appear 
nearly identical to the background, indicating little absorption during the tests. In 
Figures 40 and 41 the background scan in Figure 37 has been subtracted from the 
scans taken during the test. The variations of approximately 20 counts represent 
only approximately 1% of the signal level. In Figures 42 and 43 the background 
scan was subtracted from the average of the five scans taken at one second intervals 
and plotted in a three-dimensional representation to show the difference in trans- 
mitted intensity as a function of time and position. There was no evidence from 
the intensity scans that a laser absorption wave had been initiated, and there was 
no visual evidence from a closed circuit video camera that visible emission from a 
“water flame” had occurred. 
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IV. DISCUSSION AND RECOMMENDATIONS 


The concept of using a laser absorption wave initiated at a ceramic end-wall as 
a means to measure the laser absorption coefficient at high temperatures appears 
feasible, but was not proved in this investigation. 

Failure to initiate a laser absorption wave using the 1.5 kW laser available at 
UTSI is not surprising, since the model prediction which neglected energy losses 
due to radial conduction and convection indicated that a higher intensity would be 
needed. The experiment should be attempted again using the MSFC laser which 
should be capable of an order of magnitude greater power. Alternate methods of 
initiating the laser absorption wave might also be considered, such as interrupting 
the beam sustaining a plasma. In this case the hot decaying plasma might provide 
the initiation for a laser absorption wave. 

The use of time-of-flight measurements of the decaying blast wave from a laser 
induced spark was verified, and should provide a non-intrusive technique for mea- 
surement of translational temperature in molecular gases for which equilibrium 
between the translational and internal degrees of freedom is not a priori insured. 
Use of an excimer laser having high repetition rates to produce the spark breakdown 
insures reasonable temporal resolution for transient experiments. 

The one-dimensional model employed in this investigation should be modified 
to include two effects which are likely to be important in the initiation of a laser 
absorption wave: radial conduction in the ceramic end-wall and in the wave itself, 
and convection through the wave due to buoyancy and the expansion of the wave 
material due to heating. Both of these effects represent energy loss mechanisims, 
and would likely result in larger predicted intensity thresholds for wave initiation. 
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CONDITIONS: P = 10ATM 

20% H 2 0 
A t = .01 SECS 
l 0 = 375 W/CM 



X (CM) 

FIGURE 5 Calculated laser absorption wave for the case ‘20% HjQ in hydrogen at 10 atmo- 
spheres pressure. The incident laser intensity is 375 WJcjmr. 
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incident laser intensity is 180 W 
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AXIAL DISTANCE x (cm) 




FIGURE 8 Schematic ot the experimental arrangement for the Optoacoustic laser-beam 



- 33 - 






0. 20—i K , , . 0. 0 



- 35 - 


FIGURE 10 Calculated equilibrium concentrations for a mixture of 20% H>0 and 80% Ar 
at a total pressure of 10 atm. 



FIGURE 11 Calculated equilibrium concentrations for a mixture of 20% H 2 O and 80% 
at a total pressure of 1 atm . 








0 . 2 



02HX HX 


- 37 - 


FIGURE 12 Calculated equilibrium concentrations for a mixture of 20°% H-tO and 80 %H-> at 
a total pressure of 10 atm. 














FIGURE 17 Species concentrat ions resulting Irorn a non-equilibrium calculation for a mixture 
of 20% H 2 O and 80H He at a temperature of 2000 K and a total pressure of 10 atm. 









FIGURE 19 Species concentrations resulting from a non-equilibrium calculation for a rnixtun 
of 20% H>0 and 80% He at a temperature of 3000 K and a total pressure of 5 atm. 









FIGURE 21 Species concentrations resulting from a non-equilibrium calculation tor d mixture 
of 20% H 2 O and 80% He at a temperature of .5000 K and a total pressure of 10 a/m. 




FIGURE 22 Species concentrations resulting from a non-equilibrium calculation for a mix 
of 20 % HjO and 80% He at a temperature of 4000 K and a total pressure of 1 atm. 







FIGURE 23 Species concentrations resulting from a non-equilibrium calculation for a mixture 
of 20% HoO and 80% He at, a temperature of 4000 K and a total pressure of 1 atm. 







FIGURE 24 Species concentrations resulting from a non-equilibrium calculation for a mixture 
of 20% H 2 O and 80% He at a temperature of 4000 K and a total pressure of 5 atm. 





FIGURE 25 Species eoncentralions resulting from a non-equilibrium calculation for a mixture 
of 20% H 2 O and 80^ He at a temperature of 4000 K and a total pressure of 5 atm. 





FIGURE 26 Species concentrations resulting from a non-equilibrium calculat ion for a mixture 
of 20% H%0 and 80% He at a temperature of 4000 K and a total pressure of 10 atm. 















FIGURE 29 Species concentrations resulting from a non-equilibrium calculation for a mixture 
of 20*% #2 0 and 80% Hj at a temperature of 3000 K and a total pressure of 1 atm. 










FIGURE 31 Species concentrations resulting from a non-equilibrium calculation for a mixture 
of 20% 1 / 2 ^ and 80 °q// 2 at a temperature of .‘>000 K and a total pressure of 10 atm. 
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FIGURE 32 Species concentrations resulting from a non-equilibrium calculation for a mixture 
of 20% tjoO and 80% //•> at a temperature of 4000 K and a total pressure of I atm. 









FIGURE 34 Species concentrations resulting from a non-equilibrium calculation for a mixture 
of 20% H^O and 80 % Hi at a temperature of 4000 K and a total pressure of 10 atm . 
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FIGURE 36 Calculated sound speed based on non-equilibrium concentra tions for a mixture 
of 20% //?Q and 80% We at a temperature of 4000 K and a total pressure of 10 atrn . ~ 
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FIGURE 38 Average of 1000 scans of the array detector showing 
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FIGURE 39 



FIGURE 40 
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APPENDIX A 

FORTRAN 77 LISTING FOR LASER 
ABSORPTION WAVE MODEL 
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c *********************************************************** 

C PROGRAM TO CALCULATE UNSTEADY THERMAL WAVE IN TUBE 
C USES PATANKAR’S UNSTEADY FULLY IMPLICIT METHOD 

£ *********************************************************** 
c 
c 

DIMENSION X ( 301 ) , DELX < 301 ) , DX < 301 ) , T< 30 1 > .ABSORP (301 > ,BI(301 > 
DIMENSION CP ( 301 ) , S ( 301 ) ,TC ( 301 ), COEF (301 ) 

DIMENSION XCV( 301 ) , AKCELL ( 301 ) 

DIMENSION RHO ( 301 ) , A 1 ( 301 ) , A2 ( 301 ) ,A3(301 ) t B<301 ) 

DIMENSION ZI301 ),UD <301,2) ,LD( 301 ,2) 

C 

C X(N)-POSITION AT NTH GRID POINT (CM) 

C DELX ( N ) -X ( N + 1 ) -X ( N > 

C DX(N)-DISTANCE BETWEEN CONTROL VOLUME INTERFACES AT POINT N 
C T< N ) -TEMP . AT NTH GRID POINT (DEG. K) 

C RHO( N )-DENSITY OF MATERIAL (GM/CM~3) 

C B I ( N ) *BEAM INTENSITY AT NTH GRID LOC. (W/CM^2) 

C ABSORP ( N ) -ABSORPT I ON COEF . AT NTH GRID LOCATION 
C CP ( N ) "SPEC I F IC. HEAT AT NTH GRID LOCATION (CGS UNITS) 

C AKCELL ( N ) -THERMAL CONDUCTIVITY OF NTH CELL AT TEMP ,T(N) 

C S ( N ) “SOURCE TERM AT NTH GRID LOC. (W/CM^) 

C TC ( N ) -THERMAL COND AT C.V. INTERFACE TO RIGHT OF NTH LOC 
C COEF ( N ) - D ISCRET IZAT ION COEF. AT NTH GRID POINT. 

C 

C 

C ICOUNT-COUNTER FOR DATA SAVER (IE. SAVE EVERY ICOUNT PLOT) 

C 

c ************************************************************ 
C 

C NOTE i ALL UNITS ARE FULLY CGS UNITS!!!! mill 
C 

C ************************************************************ 
C 

c 

C TIME-TIME IN SECS! 

COMMON/SOLN/ T,BI t S 

COMMON/PROP/ TC, ABSORP, CP .AKCELL , RHO 

COMMON/GRID/ NPTS , X , DX , DELX , XCV , XL , NSURFACE 

COMMON/CLOCK/ TIME , DELT • 

COMMON/COEF / A1,A2,A3,B 

COMMON/TMP / Z , UD , LD 

C 

C 

C OPEN INPUT UNITS 
C 

OPEN(UNIT-l, NAME* ’GRID. DAT * , TYPE- ’ OLD * ) 

C 

C OPEN OUTPUT UNIT 
C 

OPEN(UNIT-2, NAME- * TOUT . DAT * , TYPE- * NEW’ ) 

C 

C READ INPUT DATA 
C 

READ(1,*) NPTS 
N I F-NPTS+ 1 
READ(1,*) NSURFACE 
READ(1,M ICOUNT 
READ(1,*) TL.TR 
C 

C TL-TEMP BC AT LEFT END , TR-TEMP . BC AT RIGHT END 
C 

READ(1,*) 10 
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c 

C 10= I NC I DENT LASER INTENSITY ON RIGHT END OF SYSTEM 
C 

BI ( NIF ) *10 
C 

READ (1,*) XL 
C 

C XL -LENGTH OF TUBE 
C 

READ ( 1 f * ) DELT , TSTOP 

C DELT -TIME STEP (SECS) 

C 

C NPTS-NO. OF GRID POINTS 

C NSURFACE - INDEX OF X LOCATION OF CERAMIC ABSORBER 
C 

CALL SETGRID(NIF) 

C 

C THIS ROUTINE SETS UP DX , DELX FROM X 
C 

CALL IC ( NPTS , T , TL , TR > 

C THIS ROUTINE SETS THE INITIAL CONDITIONS 
C 

TIME-0 

C INITIALIZE PROPERTIES AND COEFICIENTS 
CALL UPDATE (NIF) 

C 

C WRITE OUT INITIAL DISTRIBUTION 
C 

CALL OUTPUT (NPTS , TIME *T f BI ,S) 

NA-N I NT( TSTOP /DELT > 

NARR-NPTS-2 

C 

C TIME LOOP 

C 

C 

IPLOT-0 
DO 500 K- 1 , NA 
I PLOT - I PLOT + 1 
TIME-DELT*K 

CALL TDMA( A2 , A1 ,A3,B,T,Z,UD,LD, NARR , NIF ) 

IF ( I PLOT . EQ . ICOUNT) IPLOT-0 

IF( I PLOT . EQ .0 > CALL OUTPUT! NPTS , T I ME , T , B I , S ) 

CALL UPDATE (NIF) 

C 

C 

500 CONTINUE 

CLOSE ( UN I T* 1 ) 

CLOSE (UNIT-2) 

STOP 

END 
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SUBROUTINE SE'TGRIOl ISI2E) 

C 

Q ******************************************************* 

C 

C SUBROUTINE TO SET UP SPATIAL GRID FOR THE CALCULATION 
C NON-UNIFORM GRID SPACING IS ALLOWED 
C XCV-X LOC . * S OF CELL INTERFACES (CM) 

C X-X LOC. OF MAJOR GRID POINTS 
C XL -TOTAL LENGTH OF TUBE (CM) 

C 

C ******************************************************* 

C 

DIMENSION X(301 > , DX ( 301 > ,DELX < 301 ) 

DIMENSION XCV ( 301 ) 

COMMON/GRID/NPTS , X , DX * DELX , XCV f XL , NSURFACE 
DO 100 NN* 1 , NPTS+ 1 
READ! 1 ♦ * ) XCV ( NN > 

100 CONTINUE 

C 

C SET UP CONTROL VOLUMES OF "ZERO" THICKNESS AT BOUNDARIES 
C 

XCV(1>«0. 

XCV < 2 > -0 • 

XCV < NPTS + 1 > -XL 
XCV { NPTS ) -XL 
C 

C SET UP MAJOR GRID POINTS MIDWAY BETWEEN CELL INTERFACES 
C 

DO 200 NN-I f NPTS 
X(NN)-(XCV(NN) +XCV ( NN+ 1 ) ) /2 . 

200 CONTINUE 

C 

C FIND SPACING BETWEEN MAJOR GRID POINTS 
C 

DO 300 NN- 1 , NPTS- 1 
DELX(NN)-X(NN+1 )-X(NN) 

300 CONTINUE 

C 

C FIND THICKNESS OF INDIVIDUAL CONTROL VOLUMES 
C 

DO 400 NN-1 * NPTS 
DX(NN)-XCV(NN+1 ) -XCV < NN) 

400 CONTINUE 

RETURN 
END 
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SUBROUTINE UPDATE! NIF ) 


C 

C 

Q ********************************************************* 

C 

C THIS SUBROUTINE UPDATES THE VALUES OF TEMPERATURE DEPENDENT 
C PROPERTIES. CALL AT THE END OF EACH TIME STEP l ! I 11! ! ! ! !!! ! 

C 

Q *********************************************************** 

c 

DIMENSION X ( 301 ) , DELX { 301 ) , DX ( 301 > ,T<301 > ,ABSORP<301 ) 
DIMENSION B I ( 301 ) f CP ( 301 ) , S ( 301 ) , TC ( 301 ) , COEF ( 301 > 

DIMENSION AKCELL ( 301 ) ,XCV(301 ) ,RHO<301 ) 

DIMENSION A1 (301 ) , A2 ( 301 ) , A3 ( 301 ) , B ( 301 ) , Z < 301 ) , UD ( 301 ,2) 

DIMENSION LD(301,2> 

C 

C 

COMMON/SOL N/T , B I ,S 

COMMON/ PROP/ TC ,ABSORP , CP , AKCELL 

COMMON /GR I D/NPTS , X , DX , DELX , XCV , XL , NS UR FACE 

COMMON/CLOCK/TIME ,DELT 

COMMON/COEF / Al,A2 f A3,B 

COMMON/TMP/ Z,UD,LD 

C 

C 

C LOOP THRU ALL CELLS AND ASSIGN NEW PROPERTIES 

C 

C 

DO 100 NN-1,NPTS 
C 

C VARIABLE NSURFACE LOCATES THE SURFACE OF THE CERAMIC 
C 

IF(NN*GT . (NSURF ACE-1 ) ) GOTO 70 


C COME HERE IF CURRENT POINT IS IN THE CERAMIC 
C 

CALL TCSOL I D( T( NN ) , AKCELL ( NN ) ) 

CALL CPSOLID(T(NN) ,CP(NN) ) 

CALL ASOLID(T(NN) ,ABSORP(NN> ) 

CALL RHOSOL ID(T(NN),RHO(NN>) 

C 

GOTO 100 
C 

C COME HERE IF SCANNING THRU GAS 
C 

70 CALL TCGAS < T< NN > ♦ AKCELL ( NN ) ) 

CALL CP GAS (T(NN),CP(NN)) 

CALL AGAS ( T ( NN ) , ABSORP ( NN ) ) 

CALL RHOGAS (T(NN) 9 RHO( NN) ) 

C 

C 

100 CONTINUE 

C 

C NOW FIND BEAM INTENSITY PROFILE 
C 

CALL BPROFILE(NPTS,DX t ABSORP , B I ) 

C 

C CALCULATE SOURCE TERM USING NEW BEAM PROFILE AND PROPERTIES 
C 

CALL SOURCE ( NPTS+ 1 ) 
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L 

C FIND INTERFACE CONDUCTIVITIES 
C 

DO 200 NN*2 , NPTS 
T 1 "2*DELX ( NN- 1 ) 

T2*DX(NN-1 )/AKCELL(NN-l )+DX(NN)/AKCELL(NN) 

TC<NN)-T1/T2 
200 CONTINUE 

TC(1)-TC(2) 

TC(NPTS+1 >-TC<NPTS> 

C 

C UPDATE DISCRETIZATION COEFICIENTS FOR USE IN NEXT TIME STEP 
C 


DO 300 NN-2,NPTS-1 
NM-NN-1 

A1(NM)-TC(NN)/DELX(NM) 

A3(NM)-TC(NN+1 )/DELX<NN) 

AP0«RHO(NN)*CP<NN)*DX(NN>/DELT 

A2 ( NM ) *- 1 * ( AP0+A1 ( NM ) +A3 ( NM > ) 

B(NM>*-1*(AP0*T(NN)+S(NN)*DX<NN)> 

C 

C 

300 CONTINUE 

C 

B(1>«B<1)-A1<1>*T(1> 

B<NPTS-2)*B(NPTS-2) -A3 ( NPTS- 2 ) *T ( NPTS ) 

RETURN 

END 

SUBROUTINE BPROF I LE < NPTS ,DX » ABSORP * B I ) 

C 

C ***************************************************** 

c 

C THIS SUBROUTINE INTEGRATES THE BEER’S LAW EQUATION TO 
C DETERMINE THE BEAM INTENSITY PROFILE. 

C 

C ***************************************************** 
C 

DIMENSION DX ( 301 ) , ABSORP ( 301 ) , B I ( 301 ) 

C 

ETA-0. 

NS-NPTS+1 

C 

C BKNPTS+l )«INCIDENT BEAM INTENSITY (W/CM^2> 

C 

DO 100 NN«NPTS,1,-1 

ETA*ETA+ABSORP < NN ) *DX < NN ) 

I F ( ETA . GT . 10. ) GOTO 90 

BI(NN)“BI(NS)*EXP<-ETA) 

GOTO 100 

90 BHNN1-0. 

100 CONTINUE 

C 

C 

RETURN 


END 
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SUBROUTINE SOURCE(NIF) 

C 

C ********************************************************* 

C 

C THIS SUBROUTINE CALCULATES THE SOURCE TERM USING THE BEAM 
. C INTENSITY PROF I LE , TEMP ERATURE FIELD AND PROPERTIES. 

C 

C 

c ***** ****** ******* *************************************** 
C 

DIMENSION X ( 301 ) , DELX ( 301 > , DX < 301 ) ,T< 301 ) ,ABSORP < 301 ) 

DIMENSION BI(301),CP( 301 ) t S ( 301 > , TC ( 301 ) ,COEF(30I> 

DIMENSION XCV ( 301 ) , AKCELL ( 301 ) , R HOI 301 ) 

C 

C 

COMMON/SOLN/ T,BI f S 

COMMON/PROP/ TC ? ABSORP , CP ,AKCELL f RHO 

COMMON/GRID/ NPTS , X , DX , DELX , XCV , XL t NSURFACE 

COMMON/CLOCK/ TIME » DELT 

C 

C 

DO 100 NN«1,NPTS 
C 

C CHECK FOR BOUNDARY POINTS (ZERO THICKNESS CELLS) 

C 

I F ( DX ( NN ) . EQ . 0 ) GOTO 90 

CALL RAD(T( NN ) , PHI , NSURFACE , DX < NN > , NN ) 

C 

C PHI IS RADIATIVE LOSS FROM CONTROL VOLUME 
C ASSUME OPTICALLY THIN 

C 

C 

C NOW CALCULATE POWER ABSORBED FROM BEAM . 

C 

A* (BHNN+1 ) ~B I ( NN ) )/DX(NN) 

C 

C 

C S-NET SOURCE TERM 
C 

S ( NN ) *A-PH I 

GOTO 100 

90 S ( NN ) *0 . 

C 

C 

100 CONTINUE 

C 

C 

RETURN 

END 
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SUBROUTINE OUTP UT(K,TIME,T*BI ,S) 


*************************************************** 

THIS ROUTINE OUTPUTS THE SOLN (ARRAY T) ALONG WITH THE 
BEAM INTENSITY <BI(NN>>, AND THE SOURCE TERM (S(NN>) 

AT THE END OF EACH TIME STEP 

UNITS: 

T- DEG. KELVIN 
S- W/CM~3 
BI- W/CM~2 




DIMENSION T( 301 ) , B I ( 301 ) , S ( 301 ) 

DIMENSION X < 301 ) ,DX<301 ) ,DELX(301 > ,XCV<301 ) 
COMMON /GR I D/NPTS , X , DX , DELX , XCV , XL , NSURFACE 
WR ITE < 2 , 10) TIME 

10 FORMAT! IX TIME* ’.F12.6,’ SECS’) 

WR I T E < 2 * * ) NPTS 

DO 100 NN* 1 , NPTS 

BB-(BI<NN+l)+BI<NN>)/2. 

WRITE! 2, 20) X< NN ) , T( NN ) , B I (NN > , S ( NN ) 

20 FORMAT !lX,F12.6,3X f F12.6,3X,F12.6,3X,F13.4> 
100 CONTINUE 


RETURN 

END 
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SUBROUTINE IC < NPTS , T , TL , TR > 


******************************************************** 

THIS SUBROUTINE INITIALIZES THE TEMP FIELD AT TIME T=0 
THE INITIAL PROFILE IS ASSUMED TO BE A LINEAR VARIATION 
WITH X FROM T-TL TO T-TR . 

******************************************************** 

DIMENSION T ( 301 ) 

S* < TR-TL ) /FLOAT < NPTS-3 ) 

DO 100 NN*2 t NPTS- 1 
T<NN)-S*(NN-2)+TL 
100 CONTINUE 

T ( 1 >«TL 
T ( NPTS ) *TR 
RETURN 
END 
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SUBROUTINE TDMA < A , AM , AP f B , X , 2 , UO , LD , N , LENTH ) 


THIS SUBROUTINE USES A CROUT REDUCTION ALGORITHM TO SOLVE THE 
LINEAR SYSTEM OF EQUATIONS: 


A ( 1 > *X C 1 ) 

+ 

AP ( 1 ) *X ( 2 ) 

+ 

0 + 

0 ♦ . . . 

+ 

0 

-BID 

AM ( 2 ) *X < 1 ) 

+ 

A( 2 )*X ( 2 ) 

+ 

AP ( 2 ) *X ( 3 ) + 

0 * ... 

+ 

0 

« B ( 2 ) 

0 

+ 

AM(3)*X(2) 

♦ 

+ 

A( 3 ) *X ( 3 ) + 

AP { 3 ) *X ( 4 ) + ... 
♦ 

+ 

0 

« B ( 3 ) 
• 

0 

+ 

• 

• 

0 

+ 

, . . + 0 

♦ 

+ AM ( N ) *X < N- 1 ) 

+ 

A(N)*X(N) 

• 

« B ( N ) 


THREE WORKSPACE ARRAYS ARE REQUIRED UD, LD AND 2 AS WELL AS THE 
ARRAY X IN WHICH THE SOLUTION IS RETURNED. THE ARRAYS UD AND LD 
MAY ALSO BE PASSED AS X AND Z FROM THE CALLING PROGRAM, IN WHICH 
CASE THE SOLUTION OF X WILL BE PLACED IN THE FIRST COLUMN OF THE 
ARRAY UD. IF THIS IS DONE UD MUST BE PASSED AS X AND LD AS 2 OR 
THE CORRECT SOLUTION WILL NOT BE FOUND. 

NOTE: HERE N IS EQUAL TO NO. OF POINTS -2!!!! 


REAL A , AM , AP ,B,UD,LD,X,Z 
INTEGER LENTH , NPTS , I , N 

DIMENSION A< LENTH) ,AM< LENTH) ,AP( LENTH) ,B(LENTH) ,UD(LENTH,2) , 
2 LD(LENTH,2) ,X(LENTH) ,Z ( LENTH ) 

SAVE 1 -X ( 1 ) 

SAVE 2 -X ( L ENTH- 1 ) 

LD ( 1 , 1 >«A( 1 ) 

UD( 1 ,2)-AP( 1 >/LD< 1 * 1 ) 

C 

DO 50 I -2 , N- 1 
LD(I,2)-AM(I) 

LD ( I , 1 ) -A < I ) - LD ( I,2)*UD(I-1 ,2) 

UD ( I ,2 )*AP { I ) /LD ( I , 1 ) 

50 CONTINUE 
C 

LD ( N , 2 ) -AMI N ) 

L D < N , 1 )-A(N)-LD(N,2)*UD(N-l,2) 

C 

Z< 1 )»B( 1 )/LD< 1,1) 

DO 100 1-2, N 

Z( I >-<BU)-LDCI ,2)*Z( 1-1 ) )/LD< 1 ,1) 

100 CONTINUE 

C 

X ( N ) -Z ( N ) 

DO 150 I-N-1,1,-1 

X( I ) -Z < I )-UD( 1 , 2 > *X < 1 + 1 ) 

150 CONTINUE 

DO 300 J J-N , 1 , - 1 

X(00 + 1 )-X(00> 

300 CONTINUE 

X( 1 )-SAVEl 

X(LENTH-1 ) -SAVE 2 


RETURN 

END 


A-10 



on noon ooooo 


SUBROUTINE COMMENT 

Q *************************************************** 

C THIS SET OF SUBS GIVES PROPERTIES FOR 20X , 10 ATM 
C H2/H20 MIXTURES. TO ' USE : LINK AS FOLLOWS 
C SLINK TVAVE, OUTPUT, I C , TDMABV , T2010ATM 
C 

c *************************************************** 
RETURN 

END 

SUBROUTINE TCSOL ID (T ,AK) 

DATA FOR 100X DENSE POLYCRYSTALLINE AL203U! 

UNITS tW/CM-K 

DIMENSION TDATA < 6 > , PDATA ( 6 ) 

DATA TDATA/ 100. ,310.8,477,4,700. , 1033. ,8000./ 

DATA PDATA/.34, .33, .204, .130, .099, .099/ 

CALL GETPROP ( TDATA , PDATA , 6 , T , AK ) 

RETURN 
END 

SUBROUTINE CPSOLID(T,C) 

DATA FOR 100X DENSE POLYCRYSTALLINE AL203UI 
UNITS: U/G-K 

DIMENSION TDATA < 5 ) , PDATA < 5 > 

DATA TDATA/100. ,366.3,700. ,1644. ,8000./ 

DATA PDATA/. 837, 1 .15,1.42,1.51,1.51/ 

CALL GETPROP ( TDATA , PDATA , 5 , T , C ) 

RETURN 
END 

SUBROUTINE ASOLID(T,A) 

A-230. 

RETURN 
END 

SUBROUTINE RHOSOL I D ( T , R ) 

UNITS: G/CM^S 
R-3.97 

RETURN 

END 
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SUBROUTINE TCGAS<T,TC> 


C 

C DATA FOR H2 , P-10 ATM!!!!!!! 

C 

DIMENSION TDATA! 13) ,PDATA< 13) 

DATA TDATA/100. ,800. , 1600. ,2400. ,2600. ,3000. ,3400. ,3800. , 

1 4200. ,4 500. ,4600. ,4800. ,5000./ 

DATA P DATA/ 6 . 69E-04 , 3 . 35E-03 , 5 . 85E-03 , 8 . 36E-03 , 9 . 62E-03 , 

1 1 .67E-02,2.84E-02,4.6E-02,6.36E-02,7.03E-02,6.94E-02, 

2 6 . 6 1 E-02 , 5 . 94E-02/ 


CALL GETPROP ( TDATA , PDATA , 10, T, TC ) 

RETURN 

END 

SUBROUTINE CPGAS!T,C> 

C ********* ***************************** ************** 

C CP FOR H2/H20 MIX. P-10ATM. RESULTS OF NASA CODE 
C 20% H20 BY VOL 

C ************** ************** ************************* 
C“59.26*EXP(-8. 56E-07* ( T-4500 ) **2 ) +5 . 74 


RETURN 

END 

SUBROUTINE AGAS ( T , A ) 

COMMON/SP EC I ES/PTOTAL ,XH20,X0H 

q ******************************* 

C CALCULATES ABSORPTION FOR PTOTAL ATM, XH20 MOLE FRAC. H20 
C 

C SUB H20 CALCULATES MOLE FRAC. H20,0H AND SETS PRESSURE, PTOTAL, 
C 

C UNITS* CM^-1 

C DISSOCIATION IS TAKEN INTO ACCOUNT BY SUB H20. 

C 

c ********************************************************* 

CALL H20!T> 

P-PTOTAL 

ETA-XH20 

IF(T.LT. 1000) GOTO 100 
AK-8.0E-8*! T-700. >**2. 

GOTO 500 

100 AK-9 . 222E-06* ! EXP ( 6 . 66 1 5E-03*T )-l ) 

500 CONTINUE 

C 

C AK IS ESSENTIALLY ABSORP./ AMAGAT !!! 

C 

A-ETA*AK*273 . *P /T 

RETURN 

END 
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SUBROUTINE H20(T) 

C 

C ******************************************************* 

C THIS SUB CALCULATES MOLE FRACS FOR H20,OH AND SETS THE 
C TOTAL PRESSURE. MUST BE TAILOR MADE FOR EACH CONCENTRATION 
C YOU WISH TO RUN. 

C CURRENT CONFIGURATION IS 20X , 10 ATM . ! ! ! 

Q ******************************************************** 

COMMON /SP EC IES/PTOTAL , XH20 , XOH 
PTOTAL * 10 . 

C CALCULATES MOLE FRAC H20 FOR 20XH2O, P-10 ATM. 

C ALSO CALCULATES MOLE FRC FOR OH III ! 

C CURVE FIT USING DATA FROM NASA PRINTOUT. 

C 

IFIT.LE.4250. ) GOTO 100 
XH2O«.02291*EXP(-.0l694MT-4250> ) 

GOTO 200 

100 XH20-.2M2.-EXP(5.839E-26*T**6.9>) 

200 CONTINUE 

XOH«.06*EXP<-1 .28E-06*<T-4500>**2. ) 

C XOH-MOLE FRAC 0H f XH20-M0LE FRAC H201U 
C PTOTAL-TOTAL PRESSURE < ATM) 

C 

END 

SUBROUTINE RHOGAS ( T , R ) 

Q *********************************** 

C GIVES DENSITY FOR 20X H20/VOL AT 10 ATM TOTAL PRESS. 

C UN 1 TS * G/CM^3 

C CURVEFIT USING DATA FROM NASA PRINT OUT. 

C 

C **************************************************** 

c 

c 

P-10. 

c 

c HERE R IS GAS CONSTANT FOR H20 MIXTURE tATM CM^3/(GM K ) ) 

C 

I F ( T . LE . 5600 > R-18.25*EXP(-5. 73E-07* ( T-5600 ) **2 . ) + 15.75 

I F ( T . GT . 5600 ) R-34. 

C BELOW HERE R IS DENSITY <GM/CM~3> 

R-P / ( R*T ) 

RETURN 

END 
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SUBROUTINE RAD ( T , RA , NS , DXX , KK ) 

C ***************************************************** 

C THIS SUB CHOOSES THE CORRECT RAD. TERM DEPENDING ON 
C LOCATION, IE IN CERAMIC, NOT IN CERAMIC, ETC... 

C UNITS FOR RAD ARE W/CM~31!!l 

C ****************************************************** 
IF ( KK . LE . ( NS- 1 ) ) GOTO 100 

C 

C COME HERE IF IN GAS 
C 

CALL RADH20( T , RA ) 

GOTO 1000 

100 IF(KK.EQ. (NS-1 ) ) GOTO 200 

C COME HERE IF BELOW CERAMIC SURFACE 
RA-0. 

GOTO 1000 

C COME HERE IF FIRST CELL OF CERAMIC SURFACE. 

200 RA-4 .56E-12*T**4./DXX 

1000 CONTINUE 

RETURN 

END 
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ir 



SUBROUTINE GETPROP C TDATA , PDATA , NS IZE , TT , PP > 

C 

Q it * ** * ************************************************* 

c 

C INTERPOLATION ROUTINE USING LOOK UP TABLES GIVEN BY 
C TDATA * PDATA . HERE P DATA* FUNCTION { TDATA ) , 

C 

Q ***************************************************** 

C 

DIMENSION TDATACNSIZE) , PDATA < NS IZE ) 

IFCTT.LT. TDATAU)) GOTO 500 
C 

C LOOKUP TABLE LOOP 
C 

DO 100 NN«1,NSIZE~1 
NS-NN 

I F ( TT . GT . TDATA ( NN ) ) GOTO 50 
GOTO 100 

50 IF ( TT.LE. TDATA (NN+1 ) ) GOTO 200 
100 CONTINUE 

GOTO 600 

C 

C VALUE IS IN TABLE, INTERPOLATE TO GET PROPERTY 
C 

200 CALL INTER ( TDATA ( NS ) , PDATA < NS ) , TDATA ( NS + 1 ) ,PDATA(NS+1 ) ,TT,PP) 

GOTO 1000 
C 

C BACK EXTRAPOLATE 
C 

500 CALL INTER (TDATAU > , PDATA ( 1 ) , TDATA < 2 ) , PDATA ( 2 ) ,TT , PP ) 

GOTO 1000 
C 

C FORWARD EXRAPOLATE 
C 

600 Kl-NSIZE-1 

K2-NSIZE 

CALL I NTER ( TDATA C K 1 ) , P DATA < K 1 ) , TDATA ( K2 ) , P OATA ( K2 ) , TT , P P ) 

1000 CONTINUE 

C 

C NEGATIVE PROPERTIES ARE NOT ALLOWED!! 

C 

IFCPP.LT.0.) PP-l.E-12 
END 

SUBROUTINE INTER CXI,Y1,X2,Y2,X,Y) 

SL0PE*<Y2-Y1 >/(X2-Xl > 

Y*SLOPEMX-Xl )+Yl 

RETURN 

END 
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SUBROUTINE RADH20< T , RAD ) 

C 

C ********************************************************* 

C THIS SUB CALCULATES RAD . FROM H20 SPECIES GIVEN THE PRESSURE, 

C FRACTION OF H20 (AT THE GIVEN TEMP). CS(N) IS CROSS SECTION 
C CALCULATED FROM KEMP’S REPORT USING THE .1 ATMFT CURVE. 

C *********************************************************** 

c 

DIMENSION TDAT< 14) ,CSOH< 14) 

DIMENSION TDATAI 1 3 ) , CS ( 1 3 ) 

C 

C CSOH-CROSS SECTION FOR OH RAD . , CS-CROSS SECTION FOR H20 RAD. 

C 

COMMON /SPEC I ES/P TOTAL , XH20 , XOH 

DATA TDAT/298. ,500. , 1000. , 1 100. ,1500. ,2000. ,2500. ,3000, , 

1 3500. ,4000. ,4500., 5000. ,5500. ,6000./ 

DATA CSOH/8.01E-24,4.93E-23,4.93E-22,5.92E-22,6.72E-22,6.28E-22, 

1 4.26E-22,2.69E-22, 1 .88E-22, 1 . 6 1 E-22 , 1 .07E-22 , 7 . 1 7E-23 , 4 . 68E-23 , 

1 2.96E-23/ 

DATA TDATA/298 . , 500. , 1000. , 1 500. , 2000. ,2500. ,3000. ,3500. , 

1 4000. ,4500. ,5000., 5500. ,6000./ 

DATA CS/9.32E-22,2.23E-21,3.8E-21 , 3 . 35E-2 1 , 1 .97E-21 , 1 .56E-21 , 

1 1 .21E-21 ,9.38E-22,7. 1 5E-22 , 6 .03E-22 , 4 . 69E-22 , 3 . 69E-22 , 3 . 22E-22/ 

C 

C 

ANH20«XH20*PT0TAL/< 1 .362E-22*T> 

CALL GETPROP ( TDATA ,CS , 13 , T ,ACS ) 

RADC-ACS*ANH20 

RAD-2 . 268E- 1 1*RADC*T**4 . 

CALL GETPROP (TDAT,CSOH, 14, T,CST) 

ROH-l .665E1 l*PTOTAL*XOH*CST*T**3. 

RAD-RAD+ROH 

C 

C UNITS* RAD- ( W/CM^3 ) 

C 

RETURN 

END 
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SUBROUTINE TCSOL I D ( T , AK ) 


C DATA FOR 100X DENSE POLYCRYSTALLINE AL203IU 
C 

C UNITSsW/CM-K 
C 

DIMENSION TDATA ( 6 ) , P DATA! 6 ) 

DATA TDATA/100. ,310.8,477.4,700. , 1033. ,8000./ 

DATA PDATA/.34, .33, .204, .130, .099, .099/ 

CALL GETPROP ( TDATA , P DATA, 6 , T ,AK ) 

RETURN 

END 

SUBROUTINE CPSOLID(T.C) 

C 

C DATA FOR 100X DENSE POLYCRYSTALLINE AL203IM 
C UNITSt 0/G-K 
C 

DIMENSION TDATA { 5 ) , PDATAI 5 ) 

DATA TDATA/100. , 36G . 3 , 700 . , 1 644 . ,8000./ 

DATA P DATA/ .837,1.15,1.42*1.51,1.51/ 

CALL GETPROP ( TDATA , PDATA , 5 ,T , C > 

RETURN 

END 

SUBROUTINE ASOLID(T.A) 

A-230. 

RETURN 

END 

SUBROUTINE RHOSOL I D( T ,R ) 

C 

C UNITS «G/CM^3 
R-3.97 

RETURN 

END 

SUBROUTINE TCGAS(T.TC) 

C 

C DATA FOR H2, P-1 ATM! 111!!! 

C 

DIMENSION TDATA! 12), PDATAI 12) 

DATA TDATA/ 100. , 400 ., 1700 ., 2200 . ,2500. , 2800. ,3500. , 3700. , 
1 3900. ,4700. ,5000. ,14500./ 

DATA P DATA/67 .6E-05, .00121 , .0059 , .0084 , .01 3 , .022 , .068, .080 
1 , .083, .033, .024, .0565/ 

CALL GETPROP ( TDATA , P DATA , 1 2 , T , TC ) 

RETURN 

END 
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SUBROUTINE CPGAS(T.C) 


C CP FOR H2/H20 MIX. P-IATM. RESULTS OF NASA CODE 
C 20X 

C C-82.*EXP<-1 .435E-06MT-3850. )**2 >+5.704 


RETURN 

END 

SUBROUTINE AGAS(T.A) 

COMMON /SPEC I ES/PTOTAL , XH20.XOH 

C ******************************* 

C CALCULATES ABSORPTION FOR PTOTAL-1 ATM.20XH2O BY VOL.! 

C UNITSi CM~-1 

C DISSOCIATION IS TAKEN INTO ACCOUNT BY SUB H20. 

C 

Q ********************************************************* 

CALL H20< T ) 

P-PTOTAL 

ETA-XH20 

IFCT.LT. 1000) GOTO 100 
AK-8.0E-8M T-700 . )**2. 

GOTO 500 

100 AK-9 • 222E-06* < EXP < 6 * 66 1 5E-03*T ) - 1 > 

500 CONTINUE 

A«ETA*AK*273.*P/T 
RETURN 
END 

SUBROUTINE H20(T) 

COMMON/SP EC I ES/PTOTAL ,XH20,X0H 
PTOTAL-l .0 

C CALCULATES MOLE FRAC H20 FOR 20%H2O # P-1 ATM, 

C ALSO CALCULATES MOLE FRC FOR OHISU 
C CURVE FIT USING DATA FROM NASA PRINTOUT, 

C 

I F ( T , LE . 3950 . ) GOTO 100 
XH2O-.01533*EXP(- . 0207MT-3950 . ) ) 

GOTO 200 

100 XH20-,2*(2.-EXP(6. 802 E- 1 9*T* *5 * ) > 

200 CONTINUE 

XOH- .04*EXP(-2.51E-06*(T-3660>**2. ) 

END 
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SUBROUTINE RHOGAS ( T , R > 


£ ********************* ************** 

C GIVES DENSITY FOR 20X H20/VOL AT I ATM TOTAL PRESS. 

C UNITStG/CM^3 

C CURVEFIT USING DATA FROM NASA PRINT OUT. 

C 

Q **************************************************** 

c 

c 

P«1 .0 

Sr*' 

I F < T . L E . 4000 ) R»15.75*EXP( (T/4350. )**6.50> 

IFIT.GT. 4000) R«.0139*T-28.02 
R*P / ( R*T > 

RETURN 

END 

SUBROUTINE RAD < T , RA , NS , DXX , KK ) 

IF(KK.LE.(NS-1>) GOTO 100 
C 

C COME HERE IF IN GAS 
C 

CALL RADH20 ( T , RA ) 

GOTO 1000 

100 IFIKK.EQ. (NS-l ) > GOTO 200 

C COME HERE IF BELOW CERAMIC SURFACE 
RA»0. 

GOTO 1000 

C COME HERE IF FIRST CELL OF CERAMIC SURFACE. 

200 RA*4.56E-12*T**4./DXX 

1000 CONTINUE 

RETURN 

END 
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SUBROUTINE GETPROP ( TDATA , PDATA , NS I ZE , TT , P P ) 

C 

Q ****************************************************** 

C 

C INTERPOLATION ROUTINE USING LOOK UP TABLES GIVEN BY 
C TDATA, PDATA. HERE PDATA=FUNCTION(TDATA) . 

C 

Q A**************************************************** 

c 

DIMENSION TDATA < NS IZE ) , PDATA ( NS IZE ) 

IFITT.LT. TDATA ( 1 ) ) GOTO 500 
C 

C LOOKUP TABLE LOOP 
C 

DO 100 NN = 1 , NS IZE- 1 
NS-NN 

I F ( TT . GT . TDATA ( NN ) ) GOTO 50 
GOTO 100 

50 I F { TT . LE . TDATA ( NN+ l ) ) GOTO 200 
100 CONTINUE 

GOTO 600 

C 

C VALUE IS IN TABLE, INTERPOLATE TO GET PROPERTY 
C 

200 CALL INTER (TDATA! NS ) , PDATA! NS) , TDATA! NS+l) ,PDATA(NS+1 ) ,TT,PP> 

GOTO 1000 
C 

C BACK EXTRAPOLATE 
C 

500 CALL I NTER ( TDATA< 1 ) , PDATA ( 1 > , TDATA < 2 > , P DATA ( 2 > , TT , P P ) 

GOTO 1000 
C 

C FORWARD EXRAPOLATE 
C 

600 Kl-NSIZE-1 

K2-NSIZE 

CALL INTER < TDATA! K1 >, PDATA! K1 ) , TDATA! K2 > , PDATA ( K2 > , TT , PP ) 

1000 CONTINUE 

C 

C NEGATIVE PROPERTIES ARE NOT ALLOWED!! 

C 

I F < PP . LT . 0 . ) PP-1 . E- 1 2 
END 

SUBROUTINE INTER! XI , Y 1 , X2 , Y2 , X , Y > 

SLOPE*<Y2-Yl )/(X2-Xl ) 

Y e SLOP E* ( X-X 1 ) +Y1 

RETURN 

END 
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SUBROUTINE RADH20I T , RAD ) 

C 

c ********************************************************* 

C THIS SUB CALCULATES RAD. FROM H20 SPECIES GIVEN THE PRESSURE, 

C FRACTION OF H20 (AT THE GIVEN TEMP). CS(N) IS CROSS SECTION 
C CALCULATED FROM KEMP’S REPORT USING THE .1 ATMFT CURVE. 

£ *********************************************************** 

c 

DIMENSION TDAT< 14) ,CSOH< 14) 

DIMENSION TDATAt 13) ,CS< 13) 

COMMON/ SP EC IES/PTOTAL,XH20,XOH 

DATA TDAT/298. ,500. , 1000. ,1 100. , 1500. ,2000. ,2500. ,3000. , 

1 3500. ,4000. ,4500. ,5000. ,5500. ,6000./ 

DATA CSOH/8.01E-24,4.93E-23,4.93E-22,5.92E-22,6.72E-22,6.28E-22, 

1 4.26E-22,2.69E-22,1 .88E-22, 1 .6 IE-22, 1 .07E-22 , 7 . 1 7E-23 , 4 . 68E-23 , 

1 2.96E-23/ 

DATA TDATA/298 . ,500. ,1000. ,1500. ,2000. , 2500 . , 3000 . , 3500 . , 

1 4000. ,4500. ,5000. ,5500. ,6000./ 

DATA CS/9.32E-22,2.23E-21,3.8E-21 , 3 . 35E-2 1 , 1 . 97E-2 1 , 1 .56E-21 , 

1 1 .21E-21 , 9 . 38E-22 , 7 . 1 5E-22 , 6 .03E-22 , 4 . 69E-22 , 3 . 69E-22 , 3 . 22E-22/ 

ANH20»XH20*PT0TAL/( 1 .362E-22*T) 

CALL GETPROP ( TDATA ,CS, 13, T, ACS) 

RADC-ACS*ANH20 
RAD»2.268E-ll*RADC*T**4. 

CALL GETPROP ( TDAT, CSOH , 14,T,CST) 

ROH-1 .665E1 l*PTOTAL*XOH*CST*T*-3. 

RAD-RAD+ROH 

UNITS* RAD-(W/CM^3) 

RETURN 
END 


A-21 


293 


300 
20 
10 
40 
293 
500 . 
75 
. 1 


40 


0 

0 

.277778 
.555556 
.833333 
1.11111 
1 .38889 
1 .66667 

I . 94444 

2.22222 

2.5 

2.77778 

3.05556 

3.33333 

3.61111 

3.950 

4.300 

4 . 7000 

4.900 

5 

5.25 

5.5 

5.75 

6 

6.25 

6.5 

6.75 

7 

7.25 

7.5 

7.75 

8 

8.25 

8.5 

8.75 

9 

9.25 

9.5 

9.75 

10 

10.25 

10.5 

10.75 
11 

II . 25 

11.5 

11.75 
12 

12.25 

12.5 

12.75 

13 

13.25 

13.5 

13.75 

14 

14.25 

14.5 

14.75 

15 

15.25 

15.5 

15.75 

16 

16.25 

16.5 

16.75 

17 

17.25 

17.5 

17.75 

18 

18.25 

18.5 

18.75 


A-22 


c ' 3- 



19 

19.25 

19.5 

19.75 

20 

20.25 

20.5 

20.75 
21 

21.25 

21.5 

21.75 
22 

22.25 

22.5 

22.75 

23 

23.25 

23.5 

23.75 

24 

24.25 

24.5 

24.75 

25 

25.25 

25.5 

25.75 

26 

26.25 

26.5 

26.75 

27 

27.25 

27.5 

27.75 

28 

28.25 

28.5 

28.75 

29 

29.25 

29.5 

29.75 

30 

30.25 

30.5 

30.75 

31 

31.25 

31.5 

31.75 

32 

32.25 

32.5 

32.75 

33 

33.25 

33.5 

33.75 

34 

34.25 

34.5 

34.75 

35 

35.25 

35.5 

35.75 

36 

36.25 

36.5 

36.75 

37 

37.25 

37.5 

37.75 

38 

38.25 

38.5 

38.75 

39 

39.25 

39.5 

39.75 


A-23 



40 

40.25 

40.5 

40.75 

41 

41 .25 

41.5 

41.75 

42 

42.25 

42.5 

42.75 

43 

43.25 

43.5 

43.75 

44 

44.25 

44.5 

44.75 

45 


45 . 

25 

45 . 

5 

45 . 

75 

46 


46 . 

25 

46 . 

5 

46 . 

75 

47 


47 . 

25 

47 . 

5 

47 . 

75 

48 


48 . 

25 

48 . 

5 

48 . 

75 

49 


49 . 

25 

49 . 

5 

49 . 

75 

50 


50 . 

25 

50 . 

5 

50 . 

75 

51 


51 . 

25 

51 . 

5 

51 . 

.75 

52 


52 . 

.25 

52 . 

.5 

52 . 

.75 

53 


53 . 

.25 

53 . 

.5 

53 . 

.75 

54 


54 . 

.25 

54 

.5 

54 

.75 

55 


55 

.25 

55 

.5 

55 

.75 

56 


56 

.25 

56 

.5 

56 

.75 

57 


57 

.25 

57 

.5 

57 

.75 

58 


58 

.25 

58 

.5 

58 

.75 

59 


59 

.25 

59 

.5 

59 

.75 


A-24 



60 

60.25 

60.5 

60.75 
61 

61.25 

61.5 

61.75 
62 

62.25 

62.5 

62.75 

63 

63.25 

63.5 

63.75 

64 

64.25 

64.5 

64.75 

65 

65.25 

65.5 

65.75 

66 

66.25 

66.5 

66.75 

67 

67.25 

67.5 

67.75 

68 

68.25 

68.5 

68.75 

69 

69.25 

69.5 

69.75 

70 

70.25 

70.5 

70.75 

71 

71.25 

71.5 

71.75 

72 

72.25 

72.5 

72.75 

73 

73.25 

73.5 

73.75 

74 

74.25 

74.5 

74.75 

75 
75 


A-25 





UNCLASSIFIED 


5 ? AU 


REPORT DOCUMENTATION PAGE, , 

READ INSTRUCTIONS 
BEFORE COMPLETING FORM 

l. REPORT NUMBER 

2. GOVT ACCESSION NO. 

3. RECIPIENT'S CAT ALOG NUMBER 

4. TITLE (end Subtitle) 

High Temperature Measurement of Water Vapor 
Absorption 

5. TYPE OF REPORT A PERIOD COVERED 

Final Report 

6. PERFORMING ORG. REPORT NUMBER 

7. authors; 

Dennis Keefer, J.W.L. Lewis and 
Richard Eskridge 

8. CONTRACT OR GRANT NUMBERfsJ 

Contract No.: NAS8-3^320 

9. PERFORMING ORGANIZATION NAME AND ADDRESS 

The University of Tennessee Space Institute 
Tullahoma , TN 37388 

10. PROGRAM ELEMENT. PROJECT. TASK 
AREA & WORK UNIT NUMBERS 

tl. CONTROLLING OFFICE NAME AND ADDRESS 

ONRRR, Georgia Institute of Technology 
206 O'Keefe Building 
Atlanta, GA 30332 

12. REPORT DATE 

July ^11 , 1Q8S 

13. NUMBER OF PAGES 

<4. MONITORING AGENCY NAME 6 AOORESSfi/ dltteront from Controlling Ollica) 

Geo. C. Marshall Space Flight Center 
National Aeronautics & Space Admin. 

Marshall Space Flight Center, AL 35812 
Procurement Office/Vickie Carft 

15. SECURITY CLASS, (of thia report) 

UNCLASSIFIED 

15«. O EC L AS SI FlC ATION / DOWN GRADING 
SCHEOULE 

16. DISTRIBUTION STATEMENT (ot this Report) „ „ . , n n n n ), 

AP29-H* (Vicki Craft) Distribution addresses same as items 11 and lU. 

AS2UD 5 cys Cy. of transmittal letter only 

0NRR* (Michael Karp) 

EM12-3U 1 cy 
EP26 12 cys 

17. DISTRIBUTION STATEMENT (ot the abstract entered In Block 20, It different from Report) 

18. supplementary notes 


'9- KEY WOROS (Continue on reverse side It neceeemry and Identify by block number) 


20. ABSTRACT (Continue on reverse side It necessary end Identity by block number) 

* SEE ATTACHED. 

(BASA-CR-179128) HIGH lEBFEBilDfit 
P EASOREKE KT 01 kA3Efi VAPOB AESCfi ET1CN Final 
Report (Tennessee. Univ. Space Inst.) 101 p 
Avail: Nils HC AC6/BF A01 CSCL 21C 

G3/20 0082261 

. J 


N87-24525 

Unclas 


DD t^AN^ 1473 COITION OF I NOV 65 IS OBSOLETE 






